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INTRODUCTION 


This Physics Problem-Solving Set consists of a Student Lab Book and the 
corresponding Teacher’s Advisor. It was designed to help meet the need for 
computer-oriented problems in physics, providing students an opportunity to use 
a computer as a problem solving tool within a particular subject matter area. 


The materials are designed for flexible use as desired by the individual instructor. 
The material and exercises in this unit are intended as an “enrichment” ex- 
perience in the field of mechanics. 


Nearly all the exercises involve calculus in a difference approximation. This 
enables a wide variety of problems to be undertaken which would be impossible 
using only algebra and trigonometry. However, the words calculus, derivative, and 
integration are never used. The student who uses this material will have a much 
better feeling for calculus when he studies it analytically. Of course the primary 
reason a student should use the material is to better understand the ideas in 
mechanics. The use of this material will not compete with your text. Instead, it 
can be used to supplement and enrich in any fashion you choose. 


The degree of difficulty of the material is dependent upon the amount of 
assistance which you choose to provide. With no assistance, the better physics 
student should be challenged. However, given a good deal of assistance, any 
physics student should be able to work out the exercises with no great difficulty. 
The level of the material is determined by the assumption that students taking 
introductory physics will be quite capable as a group. 


The Lab Books provide text material and programming exercises for the students, 
advanced problems, and a sample program. The Teacher’s Advisor contains an 
example of a program to solve each exercise and a brief discussion of the 
important elements of exercise. 


For best results, you should study all the solutions until you are certain you have 
a complete grasp of the general methods. This should be done before assigning the 
material to the class. Generally, the exercises are cumulative so that as techniques 
are developed they are used in subsequent exercises. Therefore, you will probably 
wish to proceed through the exercises in the order in which they are given. 


The solutions should be treated as typical only. There are many ways to work a 
given problem. Students should be encouraged to develop their own solutions 
even though they may be quite different from those presented here. Also, you 
should not become obsessively concerned with computer efficiency. It is far more 
important that the student develop a logical and methodical approach to com- 
puter programming then to dwell on flashy techniques. 


MECHANICS 
Introduction 
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RATES 


Exercise 1 — Computing Velocity 


This exercise utilizes the program from Figure 1 in the Student Lab Book. The 
only change required is a new DEF statement in line 150. The intent is to let the 
student see the limiting process as the average and instantaneous velocities 
approach each other. If the student tries a value of D too small, the accuracy 
begins to deteriorate due to internal round off error by the computer. If you 
choose, you might ask the student to investigate other values of time and plot 
the results. 


100 
110 


120 
130 
140 
150 
160 
170 
180 
190 
999 


REM PROGRAM FOR AVERAGE VELOCITY 
PRINT “INPUT VALUE OF T DESIRED"; 
INPUT T 

PRINT “INPUT VALUE OF D'S 

INPUT D 

DEF FNACT)=SINCT)+Tt2 

LET V=CFNACT+D/2)-FNACT=<D/2))/D 
PRINT “AVERAGE VELOCITY IS “sv 
PRINT 

GOTO 110 

END 


RUN 


INPUT VALUE OF T DESIRED?71 
INPUT VALUE OF D?-«5 
AVERAGE VELOCITY IS 253469 


INPUT VALUE OF T DESIRED?71 
INPUT VALUE OF D?-1 
AVERAGE VELOCITY IS 24548988 


INPUT VALUE OF T DESIRED?71 


INPUT VALUE OF D?-@5 
AVERAGE VELOCITY IS 2454824 
INPUT VALUE OF T DESIRED? 1 
INPUT VALUE OF D?-@! 

AVERAGE VELOCITY IS) 245483 


INPUT VALUE OF T DESIRED?2 
INPUT VALUE OF D?el 


AVERAGE VELOCITY IS 3458463 


INPUT VALUE OF T DESIRED?2 
INPUT VALUE OF D?-@5 
AVERAGE VELOCITY IS 3458391 
INPUT VALUE OF T DESIRED?2 
INPUT VALUE OF D?-@1 


AVERAGE VELOCITY IS 3-584 
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Exercise 2 — Applications 


This exercise is based upon Exercise 1, but a very important point is developed 
here. That is, with the computer we can handle cases that would be very 
difficult even if the student had completed a course in calculus. Assuming the 
student was acquainted with differentiation he would find great difficulty with 
(c). However it is handled in the same manner as (a) and (b). 


150 DEF FNACT)= 3*T13-44*Tt2+5 
RUN 


INPUT VALUE OF T DESIRED? 1 


INPUT VALUE OF D?1 
AVERAGE VELOCITY IS) 1¢75 


INPUT VALUE OF T DESIRED?1 
INPUT VALUE OF D?e! 
AVERAGE VELOCITY [IS 1488751 


INPUT VALUE OF T DESIRED?1 
INPUT VALUE OF D?-@1 
AVERAGE VELOCITY IS) 1+86@17 


INPUT VALUE OF T DESIRED?2 
INPUT VALUE OF D?.el 
AVERAGE VELOCITY IS 286-9675 


INPUT VALUE OF T DESIRED?2 
INPUT VALUE OF D?-@1 
AVERAGE VELOCITY [S 28-8061 


156 DEF FNACT) = EXPCT)+T 
RUN 


INPUT VALUE OF T DESIRED? 1 


INPUT VALUE OF D?-1 
AVERAGE VELOCITY IS 3.71942 


INPUT VALUE OF T DESIRED?1 
INPUT VALUE OF D?-61 
AVERAGE VELOCITY [TS 3471833 


INPUT VALUE OF T DESIRED?2 
INPUT VALUE OF D?-.1l 
AVERAGE VELOCITY IS 8439216 


INPUT VALUE OF T DESIRED?2 
INPUT VALUE OF D?-@1 
AVERAGE VELOCITY IS 838928 





15@ DEF FNACT) = COSC(CTt*3)-EXPC(COS(T)) 


RUN 


INPUT VALUE OF T DESIRED?1 
INPUT VALUE OF D?! 
AVERAGE VELOCITY IS -«e633362 


INPUT VALUE OF T DESIRED?! 
INPUT VALUE OF D?.1! 
AVERAGE VELOCITY IS -1-98582 


INPUT VALUE OF T DESIRED? 1 
INPUT VALUE OF D?.-@1 
AVERAGE VELOCITY IS -1-9886@4 


INPUT VALUE OF T DESIRED?72 


INPUT VALUE OF D?e1 
AVERAGE VELOCITY IS ~-18-5488 


INPUT VALUE OF T DESIRED?2 
INPUT VALUE OF D?-61 
AVERAGE VELOCITY IS -11+2654 


INPUT VALUE OF T DESIRED72 


INPUT VALUE OF D?e@675 
AVERAGE VELOCITY TS -11.-2685 
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Exercise 3 — Automatic Accuracy 


This is an exercise in programming skill. The old value of the velocity is set equal 
to zero in line 140. The new value is computed in line 170. The test that is the 
heart of the process in line 180 gives the exit condition. Other than this, the 
program is straightforward. 


1998 


110 
126 
138 
149 
159 
166 
170 
189 


206 
210 
226 
239 
886 
881 
999 


RUN 


AT T 
AT T 
AT T 
AT T 


REM AUTOMATIC VELOCITY PROGRAM 
READ N 

DEF FNACT)=Tt2 

FOR I=1 TO N 

LET V9=6 

LET D=1 

READ T 

LET VI=CFNACT+D/2)-FNACT=-D/2))/D 
IF ABSCC(V1-V8)/V1)<-e88085 THEN 220 
LET D=D/ 19 


LET V9=V1 

GOTO 178 

PRINT “AT T = "T" V= “VI 
NEXT I 

DATA 4 

DATA 15325334 

END 

= l V = Ze 

= 2 V = Ae 

= 3 V= 6.89602 
= 4 Vs 8.-89061 
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Exercise 4 — Advanced — Plotting Results 
Hopefully, the student will discover in this exercise that for a dispiavement that is 


a linear function of time, the velocity is a constant. You might assign the students 
other linear functions and let them discover the property of constant velocity. 


120 DEF FNACT) = T 


RUN 

AT T= 1 V= 0999999 
AT T = 2 V = Le 

AT T = 3 V= le 

AT T = 4 V= le 

DONE 


Exercise 5 — Rates and Trigonometric Functions 
The student should discover that the desired function is the Cosine function. This 


is facilitated by plotting both the function for displacement and for velocity on 
the same graph. 


100 REM AUTOMATIC VELOCITY PROGRAM 


105 PRINT 
106 PRINT * T™s"V1l"*s"SINCT)%s"°COSCT)™ 
107 PRINT 


120 DEF FNACT)=SINCT) 

130 FOR T=0 TO 1 STEP «1 

140 LET vV0=0 

150 LET D=!1 

170 LET V1l=CFNACT+D/2)-FNACT=-D/2))/D 
180 IF ABSCV1-V0)<-0001 THEN 220 
190 LET D=D/10 

200 LET VO=vVI 

210 GOTO 170 

220 PRINT TeV1sSINCT) sCOSC(T) 

230 NEXT T 

999 END 
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RUN 


DONE 


V1 


i 
99498 
980079 
955343 
92113 
877619 
e 825286 
«764847 
©696659 
©621557 
540257 


~—SINCT) 


0 


9-98334E-02 


«198669 
29552 

© 389418 
© 479426 
°564643 
644218 
717356 
©783327 
841471 


COSC(T) 


le 
995004 
980067 
©955336 
921061 
877582 
825336 
«764842 
696707 
62161 
540302 








FROM RATES TO DISPLACEMENTS 


Exercise 6 — Displacement 


This exercise introduces the student to the solution of differential equations. Of 
course, he won’t know what a differential equation is; but due to the way the 


problem is phrased, the solution can be obtained without difficulty. 


100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 


999 


RUN 


REM PROGRAM TO FIND DISPLACEMENT 
REM GIVEN VELOCITY 

PRINT “INPUT INITIAL VALUE OF X “5 
INPUT X0 

PRINT “INPUT TIME INCREMENT "“S 
INPUT D 

DEF FNACT)=COSCT) 

FOR T=0 TO 1 STEP D 

PRINT TsXO 

LET X1=X0+FNACT)*D 

LET X0=X1 

NEXT T 

END 


INPUT INITIAL VALUE OF X 720 
INPUT TIME INCREMENT 72.21 


0 

el 
ee 
3 
4 
e5 
e6 
7 
8 
J 
1 


DONE 


0 

el 
©1995 
297507 
© 393041 
485147 
°572905 
©655439 
731923 
801594 
© 863755 
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Exercise 7 — Increasing the Accuracy 


The notion of a function is illustrated very nicely in this exercise. With the 
definition of the cosine function in line 160, the average velocity can be clearly 
defined in line 190. The student should be aware of the fact that different values 
of displacement are obtained using the average velocity, and that these values are 


more accurate. 


100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
999 


RUN 


REM PROGRAM TO FIND DISPLACEMENT 
REM GIVEN VELOCITY 

PRINT “INPUT INITIAL VALUE OF X "3 
INPUT X0 

PRINT “INPUT TIME INCREMENT 3 
INPUT D 

DEF FNACT)=COS(C(T) 

FOR T=0 TO 1! STEP D 

PRINT TsX0 

LET X1=X0+CCFNACT+D)+FNACT))/2)*D 
LET X0=X}1 

NEXT T 

END 


INPUT INITIAL VALUE OF X 70 
INPUT TIME INCREMENT ?-¢1 


0 

el 
°2 
e3 
o4 
5 
6 
e7 
°& 
9 
i 


DONE 


10 


0 
9e¢97TSO02E-02 
© 198504 
©295274 
389094 
479026 
564172 
°643681 
716758 
782674 
«84077 
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Exercise 8 — Finding Displacement From Acceleration 


This exercise needs to be related to the physics text to be most meaningful. The 
idea here is very powerful, and the students should work at it until the strategy is 
clear. The program is clear and needs no special explanation. 


180 


149 


RUN 


REM PROGRAM TO FIND DISPLACEMENT 
REM GIVEN ACCELERATION 
PRINT “INPUT INITIAL VALUE OF X “5S 


INPUT X@ 

PRINT “INPUT INITIAL VALUE OF V "3 
INPUT V@ 

PRINT “INPUT TIME INCREMENT "; 
INPUT D 


DEF FNACT)=T 

FOR T=6 TO 1 STEP D 
PRINT TsVG6sXO 

LET V1=V0+FNACT) *D 
LET X1=X8+V8*D 

LET vé=V1 

LET X9=X1 

NEXT T 

END 


INPUT INITIAL VALUE OF X 720 
INPUT INITIAL VALUE OF V 70 
INPUT TIME INCREMENT 7.61 


0 

ol 
e2 
3 
04 
5 
06 
e7 
oS 
9 
1 


DONE 


0 0 

0 0 

01 0 
03 001 
06 e004 
el a | 
15 02 
22] ©0035 
028 e056 
© 36 084 
© 45 12 


11 





PHYSICS 


Hewlett-Packard Computer Curriculum 


Exercise 9 — Advanced — Average Values of Velocity and Acceleration 


The only change in the program is to use average values of acceleration and 
velocity. This is done in lines 210 and 220. This is a good program to use in a 
heuristic sense. Supply the student with a wide variety of acceleration functions 
with initial conditions on velocity and position. Then the computer can be used 
to great advantage to investigate the characteristics of the motion. The student 
should be made aware that with the computer nearly any type of motion can be 
studied whereas only a few types of motion can be described if only algebra is 


used. 


100 
110 
120 
130 
140 
150 
160 
170 


190 
200 
210 
220 
230 
240 
250 
999 


RUN 


REM PROGRAM TO FIND DISPLACEMENT 
REM GIVEN ACCELERATION 

PRINT “INPUT INITIAL VALUE OF X_ "3 
INPUT X0 

PRINT “INPUT INITIAL VALUE OF V "3 
INPUT VO 

PRINT “INPUT TIME INCREMENT "5 
INPUT D 

DEF FNACT)=T 

FOR T=0 TO 1 STEP D 

PRINT TsV0sX0 

LET Vi2V0+CCFNACT+D)+FNACT))/2)*D 
LET X12=X0+¢€¢€V14V0)/2)*D 

LET VO=2V1 

LET X0=X1 

NEXT T 

END 


INPUT INITIAL VALUE OF X 70 
INPUT INITIAL VALUE OF V ?0 
INPUT TIME INCREMENT 72-1 


0 

el 
2 
o3 
04 
°5 
°6 
7 
8 
°9 
1 


DONE 


12 


0 0 

e005 -00025 
02 ©0015 
©9045 °00475 
08 e011 
e125 °92125 
18 «0365 
e245 ©05775 
32 e086 
e405 © 12225 
e5 ©1675 
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Exercise 10 — Advanced — A Generalization 


The only change required is in line 190 where the FOR statement is modified. 


190 FOR T = 0 TO 4 STEP D 
RUN 


INPUT INITIAL VALUE OF X = ?9 
INPUT INITIAL VALUE OF VO 29 
INPUT TIME INCREMENT 7.61 


0) 0 0 

el e005 00025 
2 02 e0015 
03 e945 e00475 
04 008 O11 

e5 e125 002125 
6 18 ©0365 
o7 e245 005775 
8 032 9 86 

9 «405 e12e225 
1 °5 ©1675 
lel °605 e2202T5 
12 72 0289 
1.3 e845 036725 
14 098 e4585 
15 1-125 56375 
126 1.28 6 B84 

1 7 12445 °82025 
1-8 1+62 ©9735 
1.9 1-805 1214475 
2 e Oe 1-335 
2el 22205 1254525 
202 242 127765 
203 22645 2202975 
24 2288 22306 
225 32125 2260625 
206 3238 209315 
Qe 320645 3e28275 
228 392 320661 
2.9 4.205 4.06725 
3e 4.5 4.25025 
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PEWWAWWW WWW Ww 
oe © © @ @ 
omrnrauvrbwond- 


DONE 


4-805 
5.12 

5 e445 
5-78 

62125 
648 

6-845 
722 

7 +605 
Be 


4.96775 
5 e464 
5.99225 
6 «5535 
7 014875 
7 «779 
8-44525 
91485 
9.88975 
10 +67 
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Exercise 11 — Force and Motion 


This exercise is designed to promote qualitative thinking about motion problems. 


It may help to have the student plot the results. 


800 DATA 23-432 


RUN 
Tt V 

0 -4 

1 -3.e 

Qe "2 e 

3e -le 

4 “2.83 122E-06 
DONE 


800 DATA 23-451 
130 LET F = 8 


RUN 

T V 

0 <4 

l 4e 
Le i2e 
3. CO e 
4e 286 
DONE 


x 


2 

“1-55 
-“4e1 

-5 65 

-§ 20001 


16 
9e2 
24-8 
464 


15 
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Exercise 12 — Converging on Accuracy 


As the time increment in the problem decreases the student should begin to 
suspect that the velocity and position will be more accurate. You can run this 
program with some other values for force, mass, and initial conditions to give 
more experience. 


S01 DATA «121095 


RUN 
T V x 

0 -4 2 

i -2e -j.1 

2e ~-1.84774E-06 “2 ef 

3e 2e “1.3 

4e 4e 1.59999 
DONE 


&1 DATA .055s2035 


RUN 
T V Xx 

0 -4 2 

le “2. -1-05 

2e -2e83122E-06 -2el 

3 2e “1.15001 
4.00001 3299999 1279999 
DONE 
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801 DATA 
RUN 


Q 
999999 
Qe 
3-8@8002 
4-80004 


DONE 


881 DATA 
RUN 


0 

999999 
2-80002 
3-88004 
4-00007 


DONE 


81218025 


V 


-4 

-D, 

-3.52412E-@6 
1.99999 
3.99999 


085520055 


V 


-4 

-1.99995 
4e37759E-O5 
2.88804 
4.90009 


Xx 


2 
-~1.01 
-2.02 
- 1.93001 
1.95999 


X 


2 
-1-00498 
-2-00993 
-1-01489 

1.98017 


MECHANICS 
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Exercise 13 — A Variable Force 


A simple change in the program is required. Now the force is not a constant and 
must be computed each time step. This is done in line 185 in the program. It will 
be very easy for the student to respond to this simple change. You may wish to 
point out the vast difference in the analytical techniques produced by such a 
simple change. It is precisely at points such as this that the power of the computer 


becomes visible. 


100 
110 


140 
150 
160 
170 
180 
185 
190 
200 
210 
220 
230 
240 
250 
260 
270 
800 
801! 
999 


RUN 


2 

3e 

Ae 
DONE 


18 


REM NEWTON'S SECOND LAW 
READ X0sV05M 
READ DsNoL 

PRINT 

PRINT "Ts"V"% 5 "xX" 
PRINT 

LET C=N 

FOR T=0 TO L STEP D 
LET F24*C0OS(T) 

IF C<N THEN 220 
PRINT TsV0sX0 

LET C=0 

LET V1l=VO+F*D/M 
LET X1=X0+V0*D 
LET VO=V1 

LET X0=X1 

LET C=C+]1 

NEXT T 

DATA 22-451 

DATA «131025 

END 


V X 

-4 2 

- 0544982 -2 305178 
-8-26126E-02 =e 322429 
-3-03799 -1-52959 
-6-69396 “6231051 
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Exercise 14 — Advanced — A Challenge 


This problem would be very difficult by analytical methods. However with the 
computer it merely requires several tests to locate the current value of displace- 
ment which in turn determines the appropriate value of force to use. The program 
follows without difficulty. 


100 
110 
120 
140 
150 
160 
170 
180 
181 
182 
183 
184 
185 
186 
187 
190 
200 
210 
220 
230 
240 
250 
260 
270 
800 
801 
999 


REM NEWTON'S SECOND LAW 
READ X02V05M 

READ DsNosL 

PRINT 

PRINT “T*s"VU%5*X* 
PRINT 

LET C=N 

FOR T=0 TO L STEP D 
IF ABSCX0)<S THEN 185 
IF XO >= 5 THEN 187 
LET F=+4 

GOTO 190 

LET F=0 

GOTO 190 

LET F=-4 

IF C<N THEN 220 

PRINT TsV02X0 

LET C=0 

LET V1l=VOt4F*D/M 

LET X1=X0+V0*D 


LET VO=V1 
LET X0=X1 
LET C=C+1 
NEXT T 

DATA 03501 
DATA ¢e125310 
END 


19 
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RUN 


0 

°5 

l 

1-5 

2e 

2e5 

3e 

3e5 

4e 

4e5 

Se 

505 

6 e 
6-50001 
7-090001 
7¢50001 
&-00001 
8-50001 
9.00002 
9-50002 


DONE 


V 


5 

5 

5 
3 
le 
- 999999 
“3 
-S. 
“5-8 
-5-8 
-5.8 
“5.4 
-3e4 
-1-4 


e 
nN 


wn aan on © 


eee e 
WN & WH = 


2-72 
-~2 179999 
-3-08 
-5.98 
-~8.28 
“9.58 
-9-88 
“9-18 
-7.-48 
“4-78 
-1-48 

1-82 
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Exercise 15 — Advanced — Barrier Penetration 


The only difficult point in this exercise is to insure that each step produces a 
positive velocity as long as x is less than 5. A statement must be provided to test 
for this condition. If the test is failed the object is moving in the negative 
direction and will not penetrate. For your better students, this problem is a good 
one for which to develop an automatic program to replace the trial and error 
approach. 


100 REM BARRIER SIMULATION 
110 INPUT VO 

115 LET A=vV0O 

120 LET X0=0 

130 LET V1l=V0-5-00000E-02 
140 LET X1=X0+1-00000E-02*V0 
145 IF Xi>X0 THEN 150 

146 PRINT “VELOCITY TOO LOW” 
147 GOTO 110 

150 LET VOzV1 

160 LET X0=X1 

170 IF X0<5 THEN 130 


180 PRINT 

190 PRINT AsVO 
195 PRINT 

200 GOTO 110 
999 END 


Newton’s Second Law 
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RUN 


210 

10 7-05003 
75 
VELOCITY TOO LOW 
277e5 

705 2.49996 
27 
VELOCITY TOO LOW 
277225 

7-25 1-64998 
772125 

72125 1.02498 
277-05 

7-05 «199983 
277-04 
VELOCITY TOO LOW 
2772-045 
VELOCITY TOO LOW 
2720475 

71-0475 097482 


2 
DONE 








HALF STEP METHOD 


Exercise 16 — Half Step Computation 


MECHANICS 
Half Step Method 


The half step is very easy to insert. This is done in line 175 in the program. Note 
that this does throw the velocity out of step. However the increase in accuracy for 
the position is dramatic. The values of position in the printout are exact with the 


exception of the last value. The exact value for t = 4 is 2. 


RUN 


Qe 
3e 
Ae 


DONE 


REM NEWTON'S SECOND LAW 


READ X0sV0sM 

READ DosNoL 

LET Fz=e2 

PRINT 

PRI NT eT gy" ’ eeyeee 
PRINT 

LET C=N 

LET VO=V04CF/M)*D/2 
FOR T=0 TO L STEP D 
IF C<N THEN 220 
PRINT TsV05X0 

LET C=0 

LET V1l=V0O4¢F*D/M 
LET X1=X0+V0*D 
LET VO=V1 

LET X0=X1 

LET CeC+1 

NEXT T 

DATA 22-451 

DATA e1531055 

END 


V 


“39 

-1.9 
9-99982E-02 
2el 
Ael 


Xx 


2 

mle 

=D. 

le 
1.99999 


23 





PHYSICS 


eee rere eee eee eee ee re a NE 


Hewlett-Packard Computer Curriculum 


Exercise 17 — Half Step Computation 


This exercise uses the same program developed for Exercise 16 with the data 


used in Exercise 12. 


801 DATA -1210355 
RUN 
T V 
0 “3-9 
1 “1.9 
Qe 9-99982E-02 
Je 201 
Ae 4el 
DONE 


801 DATA -05s20s5 


RUN 

T V 
0 “3-95 
le -1-95 
2e 4-99974E-02 
3e 2-05 
4.00001 4-05 

DONE 


24 


x 


2 

~le 

-D. 

“le 
1.99999 


x 


2 
“le 
-2e 
-~1-00001 
1699999 


801 DATA 
RUN 


T 


0 
999999 
Qe 
300002 
4-00004 


DONE 


801 DATA 
RUN 


t 


0 

999999 
2.00002 
3-00004 
4.00007 


DONE 


©015100s5 


V 


“3-99 

-1.99 
9-99671E-03 
2-01 
4.00999 


°005220055 


V 


“3.2995 
-1.99495 

5 -04391E-03 
200504 
4200509 


X 


2 
le 
-2e 
-~1-00001 
1.99999 


X 


2 
-.999977 
-~1.99993 
-~-999889 

2.00018 
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THE HARMONIC OSCILLATOR 


Exercise 18 — Harmonic Oscillator 





The program for this exercise poses no special difficulty. When the results are 
available, see if the students can see what kind of mathematical function would 
describe them. Generally, if the results are plotted the students will see the Sine 
and Cosine curves. This leads naturally to Exercise 19. 


100 REM NEWTON'S SECOND LAW 
110 READ XOsVOs™M 

120 READ N 

13@ LET P=3.14159 

132 LET D=P/10@ 


140 PRINT 
15@ PRINT “T's Ves XxX" 
160 PRINT 


17@ LET C=N 

18@ FOR T=@ TO 2*P STEP P/188 
185 LET F=-X@ 

190 IF C<N THEN 220 
200 PRINT TsV@sX@ 
210 LET C=@ 

220 LET V1=VO+F*D/M 
230 LET X1=X0+V8*D 
240 LET V@=V! 

25@ LET X@=X1 

260 LET C=C+1 

270 NEXT T 

800 DATA 15851 

881 DATA 18 

999 END 
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RUN 


0 

©314159 
°628318 
“942477 
12625664 
1.5708 
1.88496 
22019911 
2e91327 
2082743 
3-14159 
345575 
376991 
4.08407 
4-39823 
4.71239 
9282655 
9234071 
9 °65487 
3°96983 


DONE 
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Q 


-~ 310446 
-~ 2593443 
- -820892 
- 969876 
-1-@2497 
- 97982 1 
- «837879 
= °61214 

= °32397 

~1.08758E-@3 


3251083 
622559 
°861748 
1-81857 
1.07679 
1.82971 
880903 
*644004 
341431 


' 


l 

095579 1 
-81716 
596802 
315576 


5«30496E-04 


317691 
°607827 
°8410T2 
°993926 
1-85056 
1«@B445 
859122 
62787 

332587 


1-<67247E-03 


332687 
°637648 
882931 
1.84383 





Exercise 19 — Half Step Computation 


The half step computation is in line 175. 


100 
110 
120 
130 
132 
140 
150 
160 
170 
175 
180 
185 
190 
200 
210 
220 
230 
240 
250 
260 
270 
800 
801 
999 


REM NEWTON'S SECOND LAW 
READ X0sV0sM 

READ N 

LET P=3-14159 

LET D=P/100 

PRINT 

PRINT “Ts "VU s*X* 

PRINT 

LET Cz2N 

LET VO0=VO0-CXO0/M)*CD/2) 
FOR Tz=0 TO 2*P STEP P/100 
LET F=-X0 

IF C<N THEN 220 

PRINT TsV02X0 

LET C=0 

LET VisVO+F*D/M 

LET X1=X0+V0*D 


LET VO=V1 
LET X0=X1 
LET C=C+1 
NEXT T 
DATA 13051 
DATA 10 
END 
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RUN 


T 


0 
e314159 
628318 
942477 
1.25664 
1.5708 
1-88496 
2e1991!1 
2e51327 
2¢82743 
3e14159 
3°45575 
3°76991 
4.08407 
4.39823 
4-71239 
5°02655 
534071 
5 ¢65487 
5° 96903 


DONE 


V 


015708 
© 32546 

606279 
«830266 
974833 
1.02498 
97483 

«828331 
598928 
¢ 308357 


1-S4147E-02 


340881 
«636054 
87161 

1-02379 
1.07682 
1-02448 
°870886 
©630135 
325034 


x 


1 

950915 
«807838 
© 583908 
300341 


1-55698E-02 


© 333082 
«620989 
«850687 
999014 
1.05058 
999347 
«849342 
°©614334 
316587 


1-52419E-02 


348862 
651485 
© 893047 
1.04919 
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Exercise 20 — Changing the Mass 
This exercise should be correlated with the physics textbook for best results. 


Also, the student should plot the results to assist him to better observe the effects 
of the change in mass. 


6800 DATA 152035 


RUN 
T V x 
0 “3.-14159E-03 h 
314159 -6 05796 8E-02 ©990141 
°628318 -2127278 960732 
©942477 - «18637 0912315 
125664 - 0241903 845808 
125708 - 0292774 762486 
1 «88496 - 0337 97 66396 
2219911 -~ 37659 552142 
2251327 - 407 &58 0429213 
2282743 -~e431143 ©297575 
3e14159 - e44597 «159807 
32045575 - «452029 1 e86155SE-02 
3.76991 _ 20449185 - 123223 
408407 © 0437475 ~ 026291 
439823 ~e417112 - 2397688 
4.71239 - 388482 -.524888 
5 «02655 -e352133 ~- 641987 
5 3407 1 - «308766 - «746653 
5 65487 ~e259224 - «836799 
5 -96903 ~e2)447 - 2910618 
DCNE 
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200 DATA 12932 
RUN 


T 


0) 
314159 
0628318 
942477 
1225664 
12-5708 
1 «88496 
219911 
2051327 
2682743 
3214159 
3045575 
3076991 
4908407 
4.39823 
471239 
5 -02655 
5 -3407 1 
5 65487 
596903 


DONE 


V 


-7 -85398E-03 


-1-e53575E-02 


16383 

312508 
e446529 
559217 
©644918 
699276 
0719458 
7904303 
654381 
571979 
°460997 
326758 
17576 


© 146594 
302127 
443546 
563818 
0656912 


1 
975392 
902617 
°785023 
°628179 
439618 
022847 


5 O02351E-03 
~e219775 
~ «43486 

— 629589 
~e«7T9427 

- 920642 
-1.00229 
“1.03498 
-1-01684 
- 2948531 
- 6833161 
- 676183 
= 485121 
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Exercise 21 — Advanced — The Pendulum 
The goal in this exercise is to encourage the student to derive the equations below 
which describe the motion of a simple pendulum of length L. 

Aw/At = - (g/L)é 


Ad /At = w 


Then, the solution can be obtained from the harmonic oscillator program. 
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MORE COMPLICATED FORCES 


Exercise 22 — Damped Harmonic Oscillator 


The only change required is in the line specifying the force. 


100 
110 
120 
130 
132 
140 
150 
160 
170 
175 
180 
185 
190 
200 
210 
220 
230 
240 
250 
260 
270 
800 
801 
999 


REM NEWTON’S SECOND LAW 
READ X0sV0sM 

READ N 

LET P=3-14159 

LET D=P/100 

PRINT 

PRINT “T™%s"V% 5X" 

PRINT 

LET C2=N 

LET VO=V0-CX0O/M)*(D/2)-CV0/M) *(D/2) 
FOR T=0 TO 2*P STEP P/100 
LET F=-X0-V0 

IF C<N THEN 220 

PRINT TsV0OsX0 

LET C=90 

LET Vi=VO+F*D/M 

LET X1=X0+V0*D 


LET VO=V1 
LET X0=X1 
LET C2C+1 
NEXT T 
DATA 12051 
DATA 10 
END 


MECHANICS 
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RUN 


T 


0 

314159 
628318 
©942QA4T7 
1625664 
1.5708 
1.88496 
2019911 
2e51327 
2¢82743 
32e14159 
3°¢45575 
376991 
408407 
4.39823 
471239 
9°¢02655 
9°34071 
5 °65487 
5296903 


DONE 


V 


-8-71108E-02 
-1-84299E-02 
3¢35695E-02 
6-88752E-02 


015708 
«280771 
451363 
°53801 

555616 
521031 
451085 
°361159 
264255 
°170514 


088906 


9-60094E-02 
9-30146E-02 


082862 


6«83225E-02 
5-18057E-02 


} 
95506 
«840397 
° 684351 
511266 
© 340473 
© 185967 
5-66313E-02 
-4-31658E-02 
“2112743 
-° 154183 
-~2- 171423 
--1694e1 
-- 153463 
-e 128625 
-9-93906E-02 
-~6-94284E-02 
-e041494 
-1-74373E-02e 
1-71561E-03 





MECHANICS 
More Complicated Forces 


Exercise 23 — An Experiment 


This is very good combination laboratory and computer exercise. 


Exercise 24 — Advanced — Driven, Damped Harmonic Oscillator 


As in Exercise 22, all that is required is modification of the forces. 


100 
110 
120 
130 
132 
134 
140 
150 
160 
1706 
175 
180 
185 
190 
200 
210 
220 
230 
240 
250 


270 
800 
801 
999 


REM NEWTON’S SECOND LAW 
READ X0sV059™M 

READ N 

LET P=3.14159 

LET D=P/100 

LET T=0 

PRINT 

PRINT eT tg MU? 5 eX 

PRINT 

LET C=N 

LET VO=V0-CX0/M)*(D/2)=-C(VO0/M) *(D/2)4+(2*COS(T)/M)*(D/2) 
FOR T=0 TO 2*P STEP P/100 
LET F=-X0-V0+2*COSC(T) 

IF C<N THEN 220 

PRINT TsV05sXO 

LET C=0 

LET V1=V0+4+F*«D/M 

LET X1=X0+V0*D 


LET VO0=V1 
LET X0=X1 
LET C=C+1 
NEXT T 
DATA 1501 
DATA 10 
END 
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RUN 


T 


0 
314159 
628318 
942477 
1225664 
125708 
1.88496 
219911 
251327 
2-¢ 82743 
3°¢14159 
3245575 
3°¢76991 
4.08407 
4239823 
471239 
9°02655 
5°e34071 
5°65487 
596903 


DONE 


V 


©015708 
°272497 
386542 
338567 
13793 
--182084 
--571079 
--969553 
-1.31718 
-1-56052 
-1.-65937 
-1-59135 
“1.35415 
-- 965564 
-~2461242 
«109376 
©688574 
1.21664 
1263829 
1-9085 


l 
1-«04444 
1.15026 
1.26912 
1¢35074 
1235134 
1.24018 
1.00389 
«64809 
° 196585 


~e 311768 
-e-827912 
-1.29856 
-1-67259 
-1.-90704 
“1.97227 
-1-85547 
-1.-56242 
-1-11719 
2559911 
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ORBITAL MOTION 


Exercise 25 — Orbital Motion 


This problem should cause no particular difficulty. The symbol assignment must 
be given particular attention however. The solution given contains a half step 
computation in the velocity. 


100 REM ORBITAL PROBLEM 

110 READ X1sX2sV15V2 

lii READ DsNoL 

120 «6©LET R1=C€X1t24X2t2)t1-5 
121 LET Vl=V1-CX1*D)/ CRI *2) 
122 LET Ve=V2-(X2*"*D)/CR1*2) 
125 LET C=N 

130 FOR T=0 TO L STEP D 

135 IF C<N THEN 151 

140 PRINT TsX1sX25R!l 

141 LET C=0 

151 LET X1=X1+4V1*D 

152 LET X2@=X2+Ve*«dD 

153 LET Ri=C€X1t24xX212)91.5 
154 LET V12V1-X1*D/R1 

155 LET Ve=V2-X2*D/Ril 

156 LET C=C+! 

157 NEXT T 

800 DATA 130501 

801 DATA 1-00000E-02, 1051 





DONE 


999 END 

RUN 

0 1 0 le 
el 995004 9e-98351E-02 le 
02 «980066 © 198673 le 
°3 955336 °295525 le 
04 921061 © 389425 1-00001 
°5 877582 ° 479433 1-00001 
°6 °825335 564652 1.00001 
°7 «764841 «644228 1.00002 
799999 696706 717368 1-00002 
«899999 621609 78334 1-00003 
999999 540301 841485 1-00603 
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Exercise 26 — Experimentation 


This is strictly a discovery exercise. Students should be encouraged to try a wide 
range of initial positions and velocities and plot the resultant orbits. See if they 
can predict the orbit that will result from given initial conditions. 


Exercise 27 — A Challenge 


The initial conditions place the object at rest under the influence of a central 
force. As indicated in the printout, the object begins to acquire velocity towards 
the center of force as would be expected. The unexpected part is how can the 
object arrive at a point further out on the other side than its initial distance? The 
problem is caused by the large changes in velocity that take place as the object 
nears the center of force. These large changes invalidate the whole solution. 
Students must constantly be on the alert for this type of error. 


800 DATA 1303920 
801 DATA 2.0121951-5 





RUN 
0 1 ) le 
ol ©994992 0 985051 
2 979865 0 ©940804 
03 954303 0 86907 8 
04 ©917733 ) e7 72946 
05 e 869252 0 ©656 805 
he) «807476 0 526488 
eT ©730258 0 389429 
799999 0634073 0 ©254929 
899999 °512477 0 134593 
999999 359752 0 4e31519E-92 
Lol B8-eO3883E-N2 0 5-19491E-04 
1.3 ~4 52498 0 92-6509 
164 “6.91509 0 330 «669 
1-5 “9.30498 0 805-65 
DONE 
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Exercise 28 — Advanced — A Double Force Center 


The solution is obtained by adding an additional force to that considered 


previously. One of the fixed masses is located at (W,O) with the other at (-W,O). 


This program is fun to play with. Students can discover very interesting orbits and 


effects. 


REM DOUBLE FORCE CENTER 

READ X1sX@seV1lsV2 

READ DsNoL 

READ W 

LET R1l=C(CX1-W)te+Xer2a)11-5 

LET R@=(CX14W)te4+Xerterr1.5 

LET V1l=V1-C€X1-W)*D/CR1*2)-CX14W)*D/ CRO*2) 
LET V@=Ve-CX2*D)/CR1*2)-CXO*D)/ CREO*2]) 
LET C=N 

FOR T=0 TO L STEP D 

IF C<N THEN 210 

PRINT TsX1sX2 

LET C=0 

LET X1=X1+V1%*D 

LET X@=Xe+4+VexdD 

LET R1L=CCXL-“W)te4exXeredti1e5S 

LET R@=(CX14W)t24xX2t2)7t1-5 

LET V12V1-CX1-W)*D/R1-(€X14W)*D/R2 
LET V@=V28-CX2*D)/R1-(X2*4D)/R2 

LET C=C+i 

NEXT T 

DATA 1452020515 

DATA 1-00000E-02:;1032 

DATA 1 

END 


Orbital Motion 
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RUN 


ol 
°2 
3 
04 
5 
6 


799999 
«899999 
999999 


(XD ee eee ome me et te te ee 
eeees @# ® @ @ 
wonanrnanrkh WW = 


DONE 
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125 

1.47937 
1241946 
1.32587 
1-20676 
1.07155 
929619 
«789188 
«656693 
536439 
© 430599 
© 339437 
261709 
e 195149 
© 136995 


Be 4451 1E-02 
3°¢50322E-02 
-1-32462E-02 
-6-18257E-02 
-2e111659 
-- 163308 


0 

© 146018 
e 28435 
° 398869 
«48417 
°536144 
553918 
539354 
496264 
©429703 
¢ 345035 
«247546 
© 141963 


3°22819E-02 
-7+82628E-02 
-2 187123 
-~ 2292324 
- 2392354 
”"e 486073 
- 0572639 
-°65144 





